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Abstract 

The ANTARES neutrino telescope is a large photomultiplier array designed to de- 
tect neutrino-induced upward-going muons by their Cherenkov radiation. Under- 
standing the absorption and scattering of light in the deep Mediterranean is fun- 
damental to optimising the design and performance of the detector. This paper 
presents measurements of blue and UV light transmission at the ANTARES site 
taken between 1997 and 2000. The derived values for the scattering length and the 
angular distribution of particulate scattering were found to be highly correlated, and 
results are therefore presented in terms of an absorption length Aabs and an effective 
scattering length Ag^. The values for blue (UV) light are found to be Aabs — 60(26) 
m, Agft ~ 265(122) m, with significant (~15%) time variability. Finally, the re- 
sults of ANTARES simulations showing the effect of these water properties on the 
anticipated performance of the detector are presented. 

Key words: Neutrino telescope; Undersea Cherenkov detectors; Sea water 
properties: absorption and transmission of light. 

PACS: 07.89.+b, 29.40.Ka, 42.25.Bs, 42.68.Xy, 92.10.Bf, 92.10.Pt, 95.55.Vj 
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1 Introduction 



The Antares ^ undersea neutrino telescope [1] will use an array of photomul- 
tiplier tubes (PMT) in the deep Mediterranean Sea to detect the Cherenkov 
light emitted by muons resulting from the interaction of high energy neutrinos 
with matter. The muon track is reconstructed from the arrival time of detected 
photons. The performance of the detector is therefore critically dependent on 
the optical properties of sea water, in particular on the velocity of light and on 
the absorption and scattering cross-sections. All these parameters vary with 
the photon wavelength. The relevant spectrum spans from ultraviolet to green 
(see figure 1); the Cherenkov light spectrum varies like the photomul- 

tiplier tube quantum efficiency becomes too low to probe wavelengths longer 
than 600 nm, while the glass pressure sphere that surrounds the phototube 
absorbs the light at wavelengths shorter than 320 nm. Seasonal variations 
in sedimentation in the sea water [2] might induce variations in the optical 
parameters. 

Several measurements of sea water attenuation have been performed in the 
past. The Dumand collaboration reported an attenuation length varying with 
the light wavelength and reaching a maximum value of 60 m, with 50% ac- 
curacy, for a wavelength of 500 nm [4]. The Nestor collaboration measured a 
similar behaviour. The maximum attenuation length was found at 490 nm and 
was estimated to be 55 ± 10 m [5] . The Baikal experiment found a maximum 
absorption length of about 20 m for a wavelength of 490 nm [6]. Measurements 
performed in pure water show that the maximum of the attenuation length is 
at lower values of the wavelength in this medium (~ 400 nm) [7,8]. The max- 
imum value was measured around 90 m with 40% accuracy. More recently, 
measurements of the absorption spectrum in pure water have been reported, 
with maximum absorption lengths of 160 ± 15 m at 420 nm [9] and of 225 ± 
30 m at 417 nm [10]. 

In order to reach an optimal knowledge of the light propagation properties 
at the detector site all relevant parameters concerning photon absorption and 
scattering should be measured. These parameters, described in section 3, will 
be directly measured and continuously monitored by the Antares exper- 
iment using an instrumentation line. Until now, the adopted approach has 
been to measure these parameters with in situ autonomous devices, and these 
measurements are the subject of this paper. We present time-of-flight distribu- 
tions of photons emitted from a pulsed isotropic light source and detected by 
a PMT at different distances from the source and for two wavelengths (blue 
and UV, as indicated in figure 1). Knowledge of the time-of-fiight distribu- 
tion is essential in order to reconstruct muon tracks. While this approach is 



^ http : / / antares . in2p3 . fr 
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Fig. 1. Detection efficiency as a function of wavelength. The sohd curve shows the 
PMT quantum efficiency and the absorption by the glass of the PMT, the optical 
gel and the protective glass sphere housing the PMT [25]. The dashed and dotted 
curves are calculated with a path length in water of 5 m and 30 m respectively, 
assuming a characteristic wavelength dependence of the water absorption length 
as given in [3]. The effect of scattering is not included. The bands labelled "B" 
and "UV" indicate the wavelengths at which measurements were undertaken at the 
Antares site. 



not sufficient to fully determine the differential cross section of the photon 
scattering process, the absorption length can be measured unambiguously. A 
parameterisation which reproduces the main features of the scattering process 
can be obtained, sufficient for the needs of the tracking algorithms and of the 
detector simulation. 



2 Experimental setup and measurement procedure 



The site chosen for the deployment of the detector is southeast of Toulon 
(42°50'N 6°10'E), 40 km from shore at a depth of 2475 m. During several sea 
campaigns from 1997 to 2000 we have improved the experimental setup de- 
voted to the study of the light transmission properties, and refined the analysis 
of the data. We focus, in this section, on the final experimental configuration. 
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2.1 The mooring line 



The measuring system is mounted on an autonomous mooring line anchored 
by a sinker. The hne remains vertical through the flotation provided by syn- 
tactic buoys. After deployment, an acoustic modem ^ is used to control the 
measurement from the surface ship which stays in the vicinity of the zone 
where the line was sunk. A sketch of the mooring line including information 
on approximate heights from the sea floor is displayed in flgure 2. 

The measuring system consists of 17" pressure resistant glass spheres mounted 
on two triangular aluminum frames. A set of three mechanical cables attached 
to the vertices of the two frames deflnes their separation distance. The bot- 
tom frame supports a light source sphere which contains a set of LEDs with 
their pulsers. The top frame supports a detector sphere facing the light source 
sphere, and a service sphere (cf. flgure 5). The detector sphere houses a pho- 
tomultiplier tube, the DC-DC converter which supplies the high voltage and a 
pulse height discriminator providing a timing signal for each detected photon. 
The service sphere contains a TDC, a microprocessor which controls the mea- 
surement and records the data, and a set of lithium batteries which power the 
system. An acoustic modem remote unit is located on top of the measuring 
system. 

2.2 The light source sphere 

In order to obtain an isotropic light source for two wavelengths, 6 pairs of 
LEDs were mounted on the centres of the faces of a cubic frame 3 cm on a 
side which also supports the LED pulser boards. Each pair, which includes 
a blue and a UV LED, is covered with a single 1 cm diameter diffusing cap 
consisting of glass micro-spheres embedded in epoxy. The cube is installed at 
the centre of a 17" glass sphere whose external surface has been sand blasted to 
provide extra diffusion and to remove surface ripples or roughness which can 
destroy the homogeneity of the emitted light flux. The blue or UV emission 
colour is chosen for each new acquisition by the operator; all 6 LEDs of the 
selected colour are then flashed simultaneously. 

The spectrum of the light emitted by the LEDs was measured using a spectro- 
photometre (see flgure 3). The peak wavelength and the spectrum FWHM in 
pulsed mode operation are (375 nm, 10 nm) and (473 nm, 29 nm) respectively 
for the UV and blue LED. The time distribution of photons emitted by the 
light source sphere was measured in a dark room for both colours using the 
complete system (see figure 4) and has a FWHM of about 9 ns, with a tail 

^ ATM 845/851 from Datasonic (now Benthos), www.benthos.com 
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Fig. 2. Sketch of the mooring line used for the measurements of the water transmis- 
sion properties. The figure is not to scale. 



towards longer times. The isotropy of the source was checked to be within 
±12% by measuring the light flux for different orientations of the source sphere 
with respect to the detector sphere. 
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1600 




wavelength (nm) wavelength (nm) 

(a) UV LED (b) Blue LED 

Fig. 3. Spectral emission of the blue and UV LEDs (smoothed curves, intensity in 
arbitrary units). The light and dark areas correspond to the two extreme values 
of the range of source intensities used for the measurements, with no significant 
difference in the spectra. 

2.3 The detector sphere 

A small 1" diameter photomultiplier tube^ glued on the internal surface of 
a 17" sphere detects photons emitted by the light source. The size of the 
PMT was chosen in order to limit the counting rate from deep sea optical 
background. The PMT was then selected for its speed and low transit time 
spread. Except for the PMT window, the internal surface of the detector sphere 
is blackened in order to absorb photons outside the PMT detection solid angle. 

2.4 The service sphere 

The service sphere provides a 6 kHz trigger signal which is fed to the LED 
pulsers and through a delay to a TDC (cf. figure 5). The TDC is started by 
the delayed trigger signal and stopped by the first PMT signal above the dis- 
criminator threshold. The clock of the TDC is defined by a 40 MHz quartz 
oscillator with each 25 ns period subdivided into 32 approximately equal chan- 
nels giving an average 6t = 0.78 ns time bin. The TDC range can be adjusted 
by defining its active window; during most of the measurements it was set to 

^ Photomultipher tube 9125 SA from EMI, now ETL, 
www .electron-tub es . co . uk/splash . html 
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Fig. 4. Normalised distributions of photon arrival times measured in air with high 
statistics for the complete source sphere and with the full setup. Collimators are 
located between the source and the detector spheres. 
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Fig. 5. Sketch of the acquisition system. 

1100 channels in order to accommodate the time distributions for all of the 
source-detector distances investigated. 

The TDC linearity was studied in a dark room by recording a white noise 
spectrum with high statistics; in the absence of LED flashes, the PMT stop 
signals were provided by the background created by a controlled light leak. 
A typical white noise spectrum is displayed in figure 6, showing the non- 
linearities associated with the 32 channel subdivision pattern. A negative slope 
in the data is expected, due to the fact that the TDC is single-hit (i.e. stopped 
by the first PMT signal) and thus cannot record the arrival time of further 
photons; earlier hits are favoured over later ones. The probability pi that a hit 
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is recorded in bin i given a flat distribution of arrival time is (to first order) 
Pi = RStx {1- R ~ i? 5t X [1 - (i - 1) St] where R is the rate of 

background noise. 

CO 12000 I 1 




400 420 440 460 480 500 520 540 560 580 600 

6000 - number 



4000 




2000 



Q U I I I \ I I I I I I I \ I I I I I I I \ I J 

200 400 600 800 1000 

bin number 

Fig. 6. White noise time distribution for TDC calibration. In the enlargement, the 
non-linearities of the TDC are seen. Given their high noise level, the first (not 
shown) and last few bins were not used for the recording of physics data. 

2. 5 Experimental procedure 

The cable lengths are set on board for the desired source-detector distance 
(ranging from 15 to 44 m). Two cable lengths are typically used for each 
choice of LED voltage. 

After the line has reached the sea bed, the acoustic link is established from the 
surface ship. The light intensity is adjusted via the voltage on the LED pulsers 
to give a detection efficiency of about 1 detected photon per 100 triggers for the 
shortest source-detector distance. This ensures that the PMT is working close 
to the single photoelectron regime. The detection efficiency can artificially 
increase because of intermittent luminescence bursts. The same intensity is 
then used for the longer distance. The discriminator threshold is set at a pulse 
height value of 0.3 times the amplitude of the single-photoelectron peak, in the 
valley between the peak and the noise. For a given source-detector distance, 
several data sets corresponding typically to 5x10^ triggers each were collected 
for each of the two light colours. The overall time needed to perform the 
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measurements at one source- detector distance including the drop and recovery 
of the hne is approximately 4 hours. 

The various measurements of the water optical properties at the Antares 
site from 1997 to 2000 are summarised in table 1, in chronological order. 



Source 


(in m) 


Timp of 


Comment 


Blue 


6-27 


Dec. 1997 


Different setup (cf. text) 


Blue 


24, 44 


July 1998 


2 source intensities 


Blue 


24, 44 


March 1999 


standard 


UV 


15, 24 


July 1999 


standard 


uv 


24, 44 


Sept. 1999 


standard 


Blue, UV 


24, 44 


June 2000 


standard 


Blue, UV 


24 


June 2000 


400 m above the sea bed 



Table 1 

Data recorded for the study of the water light transmission properties. The stan- 
dard configuration is: measurement of the arrival time distribution of photons from 
a pulsed isotropic source, one source intensity (same for the two source-detector 
distances ^sd), source sphere located 100 m above the sea bed (figure 2). 

The July 1998 data taken with two source intensities II and 12 are used to 
study the systematics coming from the shape of the source time distribution, 
which exhibits a slight dependence on the LED pulser voltage (16.5 V in 
the first case, 18 V in the second). The data recorded 400 m above the sea 
bed (June 2000) were compared to the other data recorded at 100 m above 
sea bed to check the water transparency dependence over this depth range, 
corresponding to the instrumented range of the Antares detector. 

As illustrated in figure 7, all the time distributions recorded (1998 to 2000) ex- 
hibit a small tail of delayed photons. This corresponds to a small contribution 
due to scattering between the source and the detector. Figure 8 illustrates, 
for June 2000 standard blue and UV data, the shape of the photon arrival 
time distribution in air (no scattering), and in water with a source-detector 
distance rfso = 24 m or rfso = 44 m. Clearly, the width of the main peak comes 
mostly from time resolution of the setup. The FWHM of spectra recorded for 
a source-detector distance of 24 m is about 10 ns, to be compared with the 
intrinsic FWHM of 9 ns of the light source. As expected, the scattering tail 
increases with a larger separation between the source and the detector. Scat- 
tering is also seen to be more significant in UV than in blue: slightly larger 
increase of the width of the peak region and higher scattering tail, in particular 
at a distance of 44 m. 
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Fig. 7. Time distributions in UV light for the two source-detector distances (24 and 
44 m) of the first June 2000 immersion. Y-axis is proportional to the number of pho- 
tons collected. The 24 m distribution is normalised to unity. The 44 m distribution 
is normalised with respect to the one at 24 m, in addition to a (44/24)^ factor, so 
the difference between the two peaks is entirely due to the exponential attenuation 
factor. 




■20 20 40 60 80 100 120 14 

time delay (in ns) 



Fig. 8. Photon arrival times for a collimatcd air time distribution and for two distri- 
butions taken in situ, with source-detector distances of 24 m and 44 m. All distribu- 
tions are normalised to unity at the peak. Top panel: June 2000 blue data. Bottom 
panel: June 2000 UV data. The origin of the time delay axis is set to for direct 
photons. 
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The setup deployed for the first immersion (December 1997) was different 
from the one described above. An 8" photomultipher tube is located at a 
variable distance (6 to 27 m) from a coUimated continuous blue LED source 
(A = 466 nm). Only the integrated intensity was recorded. The analysis of the 
data from this immersion is described in section 4.3. All subsequent immersions 
were done with the setup described in sections 2.1 to 2.4. 



3 Simulation of the experiment 



The photon propagation (hence the time distribution of photons at a distance 
R from the source) is governed by the following inherent optical parameters: 
group velocity of light in the medium Vg, absorption length Aabs, and volume 
scattering function I3{9) = l3{9)/\ct with units m~-^ • sr~-^ (where (3{9) is the 
normalised scattering angle distribution and Agct the scattering length). The 
scattering function is roughly described by the scattering length Agct and the 
average cosine of the scattering angle distribution (or asymmetry parameter) 
(cos 6*) = 27r / 13{6) cos9 d{cos9), under the assumption of a specific shape of 
the scattering angle distribution. 



3. 1 Physics of propagation of light in sea water 



For an isotropic source of photons with intensity Jq, the intensity / detected 
at a distance R from the source by a PMT with an active area A is 



where \^^^ is the effective attenuation length, extracted from the total number 
of photons (i.e. from the integrated time distributions) recorded for two source- 
detector distances. 



An approximate degeneracy reduces the number of parameters needed to char- 
acterise the time distribution of photons at a distance R from the source. In 
particular, strong correlations can be expected if trying to extract (cos 9) and 
Agct separately, while Agf^., defined as 



~ 1- (cos^) ' 



(2) 
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describes the main part of the scattering. ^ 

We describe the scattering angle distribution following the approach of Morel 
and Loisel [12]: the scattering angle distribution is expressed as the weighted 
sum of molecular and particulate scattering. The molecular scattering is de- 
scribed by the Einstein- Smoluchowski formula for pure water, 

/5^(cos 9) = 0.06225 (1 + 0.835 cos^ 9) , (3) 

which is reminiscent of the form 

/5i^y(cos 9) = -^ {1 + cos^ 9) , (4) 



commonly called Rayleigh scattering. The 0.835 factor (rather than 1) is at- 
tributable to the anisotropy of the water molecules. The particulate scattering 
is described by the Mobley et al. [13] tabulated distribution /3p(cos 9), obtained 
by averaging the similar particulate angle distributions measured by Petzold 
in very different seas [14], at a wavelength of 514 nm. 

The total normalised scattering angle distribution is of the form: 

^(cos 9) = r//5™(cos ^) + (1 - r/)^(cos 9) , (5) 



with rj the ratio of molecular to total scattering. The average cosine of the 
total scattering angular distribution is 

(cos^) = (1 - r^) X {cos9)P = (1 - 77) x 0.924 , (6) 



^ As a general property of multiple scattering [11], the average cosine of the light 
field produced by a thin narrow parallel beam after n scattering events (cos 9)n is 
related to the average cosine for single scattering (cos 6) by the relation (cos 9)n = 
(cos . The average number of scattering events undergone by a photon reaching a 
distance R from the source is n = L{R) /Xgct where L{R) is the average path length 
of these photons. If scattering is dominantly at small angle, as in natural waters, we 
have n ~ R/Xsct- Therefore the average cosine of the light field at distance R from 
the source is: {cos9)r ~ (cos^)^/'^='=*. All combinations of Agct and (cos0) that give 
the same effective scattering length 

xeflf ^ ^sct 
'"^ -ln(cos^) 

yield the same {cos9)r. In the case where {cos 9) ~ 1, the above relation becomes 
equation 2. 
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since the average cosine of Pctzold's distribution is 0.924. In natural waters, 
T) is typically less than 0.2 [12], so (cos^) is large and the use of equation 2 is 
justified. 



3.2 The Monte Carlo simulation 

A detailed Monte Carlo simulation, which includes the geometry of the ex- 
perimental setup and the optical properties of the medium, has been used to 
analyze the experimental time distributions and extract the light transmission 
parameters at the Antares site. The parameters of the Monte Carlo are the 
absorption length Aabsj the scattering length Agct, the fraction 1] of molecular 
scattering, the source-detector distances di, the origin of time for each distri- 
bution (or the time ti at which direct photons reach the detector located at a 
distance from the source) and the collection efficiency for each distribution. 

For each photon, the distance x it will travel before being absorbed is selected 
randomly from the probability distribution proportional to exp(— x/Aabs)- The 
photon's distance to its first scattering is similarly selected according to exp(— x/ Agct)- 
If the absorption distance is shorter, the photon is propagated to its point of 
absorption and stopped. Otherwise, the type of scattering (molecular or par- 
ticulate) is selected according to their respective probabilities r] and 1 — the 
photon is propagated to its point of scattering where a new photon direction 
is sampled from the appropriate angular distribution and a new scattering 
distance is drawn. This is repeated until the total length of the photon path 
reaches its absorption distance. Time distribution histograms are filled when- 
ever the photon reaches a radial distance from the source corresponding to 
one of the possible source-detector separations. Weights are applied to take 
into account the dependence of the PMT detection efficiency on the angle 
of incidence of the photon on the photocathode [15] or to study a possible 
anisotropy of the source emission. Each Monte Carlo distribution results from 
the propagation of one million photons. 

The time distribution of the emitted light pulse is taken from the one measured 
in air for direct photons (see figure 4), its angular distribution is taken as 
isotropic, and its spectrum as monochromatic (see actual spectral width in 
figure 3) at its central wavelength. 

The obtained Monte Carlo photon spectrum corresponds to a light source with 
a vanishingly small intensity, dominated by single photon events. For realistic 

intensity conditions, since the TDC is working in the single hit mode, one 
needs to correct the spectrum for multi-photon events where the first one only 
is detected. This correction depends on the average rate of detection per light 
pulse trigger, and is calculated assuming Poisson statistics. The background 
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is removed from the data spectra (see section 4.1 on the data processing) so 
Monte Carlo spectra are generated with no noise. 



4 Data analysis 



4-1 Data processing 



For each immersion, several data sets were taken with the same configuration. 
They are all fully compatible and proved the excellent reproducibility of the 
data over a period of a few hours. They are therefore combined to reduce the 
statistical noise on the data. 

Each recorded time distribution histogram is first corrected for the non-hnearity 
of the TDC with a bin by bin division by a slope-corrected white noise time 
distribution (cf. section 2.4) as shown in figure 6. 

The optical background at this site was studied in detail [16] . It consists of a 
variable biolumincsccncc component superimposed on a constant component 
due to the radioactive decay of ^'^K. During periods without bioluminescence 
bursts, it exhibits a rate R of about 0.1 kHz/cm^, contributing as a constant 
component in the time distribution through random stop signals. Biolumi- 
nescence bursts can reach rates R up to several tens of kHz/cm^ generally 
lasting for hundreds of micro-seconds to seconds i.e. for longer than the time 
range of the TDC. Given the highest rates observed, bioluminescence bursts 
only contribute as an additional noise which appears as a linearly decreasing 
component. Extrapolating the result of section 2.4, the total number of back- 
ground events Ni in bin i of the spectrum in the region free of LED events is 
therefore given by 



-^triggers 

Ni= J2 Rj St[l- {i-l)Rj 5t] 

-^triggers \ /-^triggers 

E Rjst]-{i-i)[ y: (Rj^ty 

i=i / V j=i 

= a-(i-l)6, (7) 



where A^triggers is the number of triggers and Rj is the background rate (from 
bioluminescence and the decay of ^°K) during cycle j on the 1" PMT {R = 
RS). The background contribution is determined by a first-order polynomial 
fit to the data in the region free of LED hits (i.e. before the signal from the 
LED direct photons): — a — {i — l)b. Over 500 bins are available for the noise 
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fit. The background rates were seen to vary between 0.2 and 2.5 kHz/cm^. 



When taking into account the fact that the multi-hit correction apphes to 
the sum of the hits from the background noise and from the LED source, 
the background to be subtracted from each bin i of the data is given by the 
following equation (which results from the difference between the expected 
distributions in the presence and in the absence of noise) : 



-^triggers 

-^noise(0 ^ ] 



i-l 



{Rj 6t + i?LED, 6t)[l-{i- l)Rj 6t-J2 ^LED, St 

\ fe=l 

—-^triggers X RhEBi ~ ^ -^LED^ Stj 

a - (i - 1)6 - a ^ -Rled^ St - {i - 1) a RhEDi St (8) 



k=l 



where a and b are the same as in equation 7 and i^LEDj is the rate of hits 
in bin i coming from the LED source (independent of the cycle number). 
Three correcting terms to the canonical value a of the background appear in 
equation 8. The slope b in the noise remained small {b < 10^'^). The intensity of 
the source was chosen so as to minimise the multi-hit correction (and maintain 
a reasonable signal-to- noise ratio), so that Z)fel\(i?LEDfe St) is at most ~ 5% 
and the second correcting term also remains small. The last term, however, 
can be quite large in bins where the rate of hits from the LED source is large. 

The result of this procedure is illustrated in figure 9 with the example of blue 
data taken in June 2000. The spectrum results from a total of 6 x 10^ triggers. 
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Fig. 9. Noise counts subtracted from the blue data of June 2000. The thin soUd fine 
is the level the noise would have had in the absence of LED hits. 

with an average collection efficiency of 6%. The noise level estimated from 
the first 700 bins is a = 6.05 with a slope b = -8.3 x 10"^. Bin number 750 
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contains the highest signal, with 3335 counts and a noise contribution of —9 
hits, meaning that the photons collected earlier in the spectrum prevented 9 
signal counts from reaching this bin. The noise level in the tail of the spectrum 
is stable at a level of 5.6 counts. 

After the subtraction of the background, when two different source-detector 
distances are available for the same date and source intensity, a measure of the 
effective attenuation length is obtained from the integrated time distributions 
(cf. section 4.3). The experimental time distributions are then fitted simulta- 
neously to Monte Carlo distributions to extract the absorption and scattering 
parameters (cf. section 4.4). 

4-2 Group velocity of light 

The group velocity of light can be computed from the equation 



using various empirical models for the index of refraction n evaluated for the 
parameters of the Antares site (pressure p — 230 atm, salinity S — 38.44 °/oo 
and temperature T — 13.2°C). 

Using four different experimental data sets of pure water and sea water under 
various pressures, Millard and Seaver [17] (thereafter referred to as MS) have 
developed a 27-term algorithm that gives the index of refraction to part-per- 
million accuracy over most of the oceanographic parameter range (salinity 
S = — 40 °/oo, temperature T = — 30°C and pressure p = 1 — 1080 atm), 
but only over a hmited range of wavelengths (500 — 700 nm) so that we need 
to extrapolate to use it for our wavelengths of 375 and 473 nm. The result is 
illustrated in figure 10, curve labelled MS. 

A simple empirical equation for the index of refraction of sea water n(A, S, T) 
for A in [400 — 700] nm can also be found in [18] by Quan and Pry (thereafter 
referred to as QF), based on data from Austin and Halikas [19]. The pres- 
sure dependence was not included in their equation, so we added it assuming 
the same linear dependence as that observed on pressure-temperature plots 
from [20] . The wavelength dependence of this model is the curve labelled QF 
in figure 10, showing excellent compatibility with the MS model. 

As an experimental verification, the following consistency check was per- 
formed. From the time distributions recorded with two source-detector dis- 
tances for a unique source intensity, one can extract the group velocity of 
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light Vg = Ad/ At where Ad is the difference between the source-detector dis- 
tances for the two immersions, and At is the difference between the times at 
which the non-scattered photons emitted by the source reach the detectors 
located at the two distances from the source. 

A length difference Ad = 20.60 ± 0.f4 m was measured on shore imder a 
tension of 50 kg on the cables to simulate the buoyancy pull of the immersed 
line. ^ The time differences are 



At 



94.3 ± 0.1 (stat.) ± 0.1 (syst.) ns (Blue) 
95.7 ± 0.1 (stat.) ± 0.1 (syst.) ns (UV) 



(10) 



where the first error is the statistical error and the second a systematic error 
coming from the measurement of the electrical length of the electric cables 
joining the source to the detector (a cable was associated with each source- 
detector distance). These data therefore imply the following velocities of light 
(the error includes the uncertainties on Ad and At stated above), also plotted 
in figure 10: 



Vg (experimental) = < 



0.2185 ± 0.0015 m/ns (Blue) 
0.2153 ± 0.0015 m/ns (UV) 




QF + pressure correction 
MS 



0.212 



350 400 450 500 550 600 

wavelength (nm) 

Fig. 10. Comparison of measurements of the group velocity of light with model 
predictions, for A = 472.5 nm (Blue) and A = 374.5 nm (UV). The phase velocity 
function of wavelength is also shown. 

As can be seen from the figure, experimental and analytical values are in good 
agreement. Given the large uncertainty on the determination from the data of 



^ Parafil cables were seen to stretch by about 1% under a tension of 50 kg, although 
specifications mentioned a stretch of at most 1.6 %o- Part of the uncertainty on Ad 
comes from the uncertainty on the stretchability. 
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the group velocity of light, however, the value of Vg in the Monte Carlo is set 
to the average of the analytical estimates described above, i.e.: 



It should be noted that the results of the Monte Carlo fit are not influenced 
by a small change in the chosen value of the group velocity. 

4-3 Effective attenuation length 

The effective attenuation length A^f^ gives an indication of the fraction of the 
photons emitted by the source that are detected (inchiding those that reach 
the detector although they have scattered on their way). It varies with the 
angular distribution of the source emission and with the angular acceptance 
of the detector. It can be computed from the ratio of the total light detected 
at the two source-detector distances di and 



where N^^ (t) is the time distribution at distance di after background subtrac- 
tion and multi-photon event correction. The effective attenuation lengths and 
the corresponding statistical errors for the data listed in table 1 are given in 
tables 2 and 3. 

An uncertainty in the noise subtraction procedure would not affect this result, 
as explained in section 4.5.2. The two data sets taken in July 1998 with differ- 
ent LED intensities yield compatible values of the effective attenuation lengths 
(62.6 ± 1.0 m and 60.3 ± 0.4 m), despite a large correction for multi-photon 
events in the second case. 

As stated in section 2.5, the effective attenuation length was also measured 
for the different setup deployed in December 1997, which used a continuous 
collimated source. While the distance D between the source and the PMT was 
varied from 6 to 27 m, the intensity of the source $led was adjusted so as 
to yield a constant current /pmt on the PMT. The setup was cahbrated with 
a similar procedure in air. The emitted and detected intensities in water are 
related by 




(12) 




(13) 




(14) 
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making it possible to estimate the effective attenuation length from the de- 
pendence of the required LED intensity with the distance (cf. figure 11). The 
agreement of the data with a decrease as given in equation 14 yields an effec- 
tive attenuation length 

A°ft (Blue, coUimated) = 41 ± 1 (stat.) ± 1 (syst.) m . (15) 

The coUimation of the source prevents a direct comparison with values given in 
table 2. A Monte Carlo simulation describing the two setups shows, however, 
that the above \'^^ would yield A^f^ = 44 ± 1 (stat.) ± 1 (syst.) m with the 
present (isotropic) setup, similar to the result found in June 2000. 

Q 

LU 

9 

CM 

Q 

8 
7 
6 



5 10 15 20 25 30 

Distance D (in m) 

Fig. 11. Determination of the effective attenuation length from the setup immersed 
in December 1997 (see equation 14). 

4.4 Absorption and scattering lengths 

The uncertainty on the exact cable lengths strongly affects the time of arrival 
of direct photons, but has a negligible effect both on the shape and on the 
amplitude of the time distributions. The source-detector distances di and d2 
are therefore fixed to the values measured on shore under a tension of 50 kg, 
while the direct photons arrival times ti and ^2 are unconstrained to absorb 
the uncertainties on the distances and avoid biasing the results. The group 
velocity of light is taken from equation 12. The other free parameters of the 
fit are the absorption length Aabs, the fraction of molecular to total scattering t] 
(or equivalently {cos 9) as explained in equation 6) and the effective scattering 
length Aft. 
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All the results are summarised in tables 2 and 3, and plotted in figure 12. 
The parameters Aabsj ^tct V result from the global fit described above. 



The efi^ective attenuation length A^^t is computed according to the method 
described in the previous section. Therefore, the usual equality between 1/A 



cff 
att 



and the sum of the inverses of Aabs and A|"^ does not hold, as they were derived 
from different methods. The scattering length Agct is obtained from Ag^ and rj 
according to equations 2 and 6. 

The March 99 data were recorded with a low source intensity, and thus a 
low collection efficiency, resulting in a high statistical error on the measured 
parameters. The data sets with two different light intensities available for the 
July 1998 immersion yield fully consistent results (to within la) for the fit 
parameters. Only the mean values are therefore reported here. 



Epoch 


\efr 

(in m) 


Aabs Ag^ 

V 

(in m) (in m) 


Asct 

(in m) 


July 1998 
March 1999 
June 2000 


60.6 ± 0.4 
51.9 ±0.7 
46.4 ± 1.9 


68.6 ±1.3 265 ±4 0.17 ±0.02 

61.2 ±0.7 228 ±11 0.19 ±0.05 

49.3 ±0.3 301 ±3 0.05 ± 0.02 


62 ±6 
58 ±18 

38 ±8 



Table 2 

Summary of the results for the blue data (statistical error only). 



Epoch 


\eff 
'^att 

(in m) 


\ \eff 
^abs ^sct 

V 

(in m) (in m) 


Asct 

(in m) 


July 1999 
Sept. 1999 
June 2000 


21.9 ±0.8 
22.8 ±0.3 
26.0 ±0.5 


23.5 ±0.1 119 ±2 0.16 ±0.03 

25.6 ±0.2 113 ±3 0.18 ±0.01 
28.9 ±0.1 133 ±3 0.12 ±0.01 


27 ±4 

28 ±1 
24 ±1 



Tabic 3 

Summary of the results for the UV data (statistical error only). 

All fit x^'s per degree of freedom are about 1 (within ±0.5). Figure 13 illus- 
trates, for June 2000 blue and September 1999 UV data, the photon arrival 
time distribution in water with a source-detector distance dsn = 24 m or 
dsB — 44 m, with the Monte Carlo fit superimposed on top of the two in-situ 
distributions. 



To illustrate the impact of the water transparency properties at various epochs, 
time distributions have been generated with each set of best-fit parameters 
assuming a unique setup, and in particular, a unique choice of the source time 
distribution: that of the June 2000 UV data. The distributions obtained are 
shown in figure 14. 

It can be seen in the figure that the absorption length in UV is smaller than in 
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Fig. 12. Absorption (dots) and effective scattering (triangles) lengtiis measured at 
the Antares site at various epoclis for UV and blue data. Horizontal error bars 
illustrate the source spectral resolution (it la). The large circles are estimates of the 
absorption and scattering lengths in pure sea water (from [20]). The dashed curve 
is the scattering length for pure water [21], upper limit on the effective scattering 
length in sea water. 

blue (lower relative height, in UV, of the distribution at the largest distance 
compared to that at the shortest). The higher tail of delayed photons for the 
UV distribution is compatible with a smaller scattering length in UV. The 
UV distributions show little dispersion. The first two blue distributions (July 
1998 and March 1999) are quite similar. Only the June 2000 blue distribution 
differs significantly from the other two, with a much lower tail of delayed 
photons due to its much larger effective scattering length. Its statistical error 
bars, however, are also very large, due in particular to the long tail of the time 
distribution of the June 2000 blue source in the calibration spectrum, which 
reduces the significance of the tail of delayed photons in the data. 

Effective attenuation length results from a combination of absorption and 
scattering. Since the scattering tail is small (confirmed by the large value of 
compared to Aabs), the absorption length should be close to the effective 
attenuation length. As expected, ^ Aabs in ah data sets. The correlated 
variations of A^f^ and Aabs — quantities determined independently and quite 
robustly — strengthen the hypothesis of fluctuations of the medium optical 
properties with time. 

As explained in section 3.1, A^f^. is used instead of the more physical Agct in 
order to avoid large (~ 80%) correlations between Agct and r]. The correlation 
coefficient between A^J^ and r] is typically less than ~ 10%. All other corre- 
lation coefficients are compatible with zero, confirming the independence of 
absorption and scattering properties. 

Despite the comments of section 3.1, one can try to extract simultaneously 
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Fig. 13. Distributions of photon arrival times with the best-fit Monte Carlo curves 
superimposed on top of the data. The same normalisation procedure as in figure 7 
is applied. Top panel: blue data recorded in June 2000. Bottom panel: UV data 
recorded in September 1999. 

Available separate figure : scat_art.jpg 

Fig. 14. Simulated distributions of photon arrival times assuming source-detector 
distances of 24 and 44 m and the experimental setup of June 2000 UV data , for 
the best fit parameters of each set of data. The same normalisation procedure as in 
figure 7 is applied. Top panel: Blue results, bottom panel, UV results. 

the molecular scattering length A^^ = Xsct/V: the particulate scattering length 
Agct = Asct/(1 — v): cind the shape of the particulate phase function using a 
generic one-parameter Henyey-Greenstein function (3^^[g,cosd) : 

^o^g^cose) = ^ , (16) 

An [I + — 2g cos (J y^'^ 



where g = {cos 6). The results are illustrated with the example of the Septem- 
ber 1999 UV data. As expected, a fit with these parameters yields very large 
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dot position 


star position 


Asct = 30 m 


Asct = 76 m 


J] = 0.25 


7] = 0.56 


(cos) p = 0.924 


(cos) p = 0.54 



Table 4 

Scattering parameters at two locations of the degeneracy valley indicated with a 
dot and a star respectively in figure 15. 

correlation coefficients between the various scattering parameters : 

Aabs A- A^t (cos)f 
Aabs 1 -0.19 -0.12 -0.08 

1 -0.87 -0.96 . (17) 
1 -0.84 

(cos)P 1 



\m 
^sct 

^sct 



The plots of figure 15 show the existence of a minimum as a function of 
Ag^, here found for A^^,. ~ 100 m, but indicate large degeneracy between rj 
and (cos 6*) (region delimited by the dashed line). For instance, the dot and 
the star shown in these plots have radically difTerent scattering parameters (cf. 
table 4), although the time distributions for either set of parameters are almost 
indistinguishable. This follows from identical scattering angle distributions 
except for small angles (less than ~30 degrees) to which the experiment is 
poorly sensitive. 

Available as a separate figure : EtaMcEff_NB.jpg 

Fig. 15. r] vs. (cos) for various slices in Af^. The various shades of grey separate 
equidistant regions with a difference of 0.5 in the normalised P^r degrees of 
freedom. The dashed line delimits the la valley. The dot and the star correspond 
to the scattering parameters in table 4. 

This degeneracy in the description of the scattering properties explains the use 
of the widely cited scattering measurements of Petzold [14] for the particulate 
scattering phase function, only allowing 77, the fraction of molecular to total 
scattering, to vary. 



4-5 Discussion on systematic uncertainties 



Several sources of systematic uncertainties may affect the results of the anal- 
ysis. These include the slight anisotropy of the source (see section 2.2), an 
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uncertainty in the noise determination (see section 4.1), the reproducibihty of 
the source intensity at a given voltage for the simuhaneoTis fit of time distri- 
butions recorded several hours apart (for the two source-detector distances), 
the knowledge of the time resolution of the source (time distribution in the 
absence of scattering, measured in the lab, cf. section 2.2) and the angular 
detection efficiency of the detector module. The systematics that affect the 
relative normalization of the spectra taken with the two source-detector dis- 
tances will mostly have an impact on the efi^ective attenuation length and on 
the absorption length, while those that change the shape of the spectra will 
rather affect the scattering length. These systematics are studied with the help 
of the Monte Carlo described in section 3.2. An estimate of their effect is given 
below. 



4-5.1 Source anisotropy 

To test the impact of the 12% anisotropy of the source, time distributions 
are generated with various angular distributions of the light source emission. 
While the tail of the time distribution is mostly produced by photons having 
scattered at large angles, the peak (containing most of the hits) comes from 
photons that have either reached the detector directly or scattered at small 
angles, i.e. in either case comes from photons emitted by the portion of the 
source facing the detector. Since the relevant information for the determination 
of Aabs and \'^^ is the relative normalization of the time spectra collected with 
two source-detector distances, which to first order is the relative peak height, 
and since the source sphere remains in the same position with respect to 
the detector for all the measurements, the source anisotropy has a negligible 
impact on the absorption and the efi'ective attenuation lengths. It induces a 
small distortion on the tail of the time distributions, however, with an impact 
of the order of 4% on Agct or Ag^^.. 

4-5.2 Noise subtraction 

An uncertainty in the noise subtraction will not affect the absorption nor the 
effective attenuation lengths since the noise remains small compared to the 
signal, and absorption dominates largely over scattering. On the other hand, 
the scattering length is mostly determined from the tail of the photon arrival 
time distribution, which is barely above noise level. The scattering length is 
therefore crucially dependent on the precision of the noise subtraction. For 
data taken after June 1999, the spectra have 1100 time bins. The noise contri- 
bution before the main signal peak is obtained by a fit over more than 500 bins 
with negligible impact of the uncertainty on the determination of the noise 
(< 1%). The acquisitions of July 1998 and March 1999, on the other hand, 
were done with 128 bins only. The noise is therefore estimated from only ~ 20 
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bins in the worst case (data with the short source-detector distance). The un- 
certainty on the noise determination can then induce an uncertainty on the 
effective scattering length of at most 8%: 16 m, 13 m and 3 m for the July 
1998, March 1999 and June 2000 blue results respectively, and 1 m for all the 
results in UV. 



4-5.3 Stability of LED intensity and PMT efficiency 

The LED intensity and the PMT efficiency are stable over a given immersion, 
as was verified by the excellent reproducibility of time spectra taken under 
the same conditions up to two hours apart. On two consecutive immersions 
with different source-detector distances, however, it is not possible to check 
the stability of the setup. This affects the global normalization of the time 
spectrum, i.e. the photon collection efficiency (the larger the intensity, the 
higher the collection efficiency), and thus the effective attenuation length (see 
section 4.3). Depending on the intensity actually used during each immersion, 
an assumed uncertainty of 1% on the LED intensity and PMT efficiency affects 
the effective attenuation length (and similarly the absorption length) by 1 to 
11%: 4.5 m, 5.5 m, 1 m and 2 m for the July 1998 (intensity II), July 1998 
(intensity 12), March 1999 and June 2000 blue data respectively, 2 m for the 
July 1999 and September 1999 UV data and 1 m for the June 2000 UV data. 
The impact on the scattering length is negligible since the shape of the time 
spectrum is unchanged. 

4.5.4 Source time distribution 

The shape of the time distribution of the source varies with the LED input 
voltage, as illustrated in figure 16. Except for the data recorded in June 2000, 
time spectra of the source were not always available at the exact same voltage 
as the one used for the water measurements. 

Because the relative normalization of spectra at different source-detector dis- 
tances is independent of the choice of the source intensity, the effective at- 
tenuation length, and therefore to first approximation the absorption length, 
are not affected. The convolution by a slightly different shape of the source 
spectrum does not affect the level of the scattering tail which is set by the 
effective scattering length. The latter is therefore also little affected by a poor 
knowledge of the source time spectrum. The strongest effect appears in the 
peak of the distribution, and reflects on the angular dependence of scattering, 
i.e. on the ratio of molecular to particulate scattering. More quantitatively, 
the data from July 1999 (taken with a DAC voltage of 5.15 V) are fitted with 
two different source time spectra, first with one at 4.75 V then with one at 
5.5 V. The absorption length is unchanged, the effective scattering length in- 
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Fig. 16. Normalised time spectra of the UV source for DAC voltages of 4.75V (solid 
curve), 5V (dashed curve) and 5.5V (dotted curve). 

creases by 6 m between the two fits and rj increases by 0.06. Extrapolating 
these results to the difference between the voltage used under water and the 
one used for the measurement of the source spectrum yields an uncertainty of 
3 m for the July 1998 and March 1999 bhic data and of 1.2 m and 2.4 m for 
the July 1999 and September 1999 UV data respectively There is no effect 
on the June 2000 blue or UV results since the same source voltage is used for 
calibration and data taking. 



4-5.5 Angular acceptance of the detector sphere 

The angular detection efficiency (studied in [22]) of the glass sphere housing 
the photomultiplier tube was simulated under the assumption of the exact 
knowledge of various factors, some of which are affected by large uncertainties 
as for instance the thickness of the photocathodc, the thickness of the opti- 
cal gel which ensures the optical contact between the photocathodc and the 
glass sphere, or the complex refractive index of the photocathodc. The largest 
changes in the shape of the efficiency curve are obtained by considering the 
smallest (respectively largest) value of the complex refractive index within its 
possible range (1.10 + 1.70i to 2.75 + 2.50i), together with the smallest (resp. 
largest) value of the photocathodc thickness (between 16.4 nm and 26.5 nm). 
Both of these configurations are shown in figure 17. 
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Fig. 17. Extreme models of the detector sphere angular efficiency (arbitrary units). 
The sohd curve uses a parameterisation of the Moorhead and Tanner data [15] and 
is the model used in the simulation of the experimental setup. Model a (resp. model 
b) uses the thinnest (resp. the thickest) of the photocathodes together with the 
smallest (resp. the largest) complex refractive index. 



From figure 17 it can be seen that a change in the parameters mentioned above 
reflects upon the ratio n<45o /r;,>45o of photons detected at angles smaller or 
larger than about 45°, where the efficiency shows a local enhancement. With 
model a on the one hand, this ratio decreases, causing more photons to be 
detected at large angles and therefore requiring a larger scattering length to 
reproduce the light transmission data. Model b on the other hand requires a 
shorter scattering length. In both cases, the effect is of the order of 8% on the 
effective scattering length. 



4-5.6 Summary of systematic uncertainties 

The total effect of the systematic uncertainties mentioned above is summarised 
in tables 5 and 6, adding all the systematic errors in quadrature. 

The systematic error is significantly larger than the statistical error. It results 
in an uncertainty of 5 to 11% on the light transmission parameters. Given 
these uncertainties, the variations with time of the values of the parameters 
(cf. section 4.3) are reduced to a ~ 2cr effect. To cope with possible temporal 
variations of the light transmission parameters, which in this paper are shown 
to be small, the ANTARES detector will monitor continuously the optical 
properties at the ANTARES site. 



Simulation model 




60 70 80 90 
zenithal angle (°) 



30 



H/pocn 




July 1998 
March 1999 
June 2000 


60.6 ±0.4 ±5 68.6 ±1.3 ±5 265 ±4 ±28 
51.9 ± 0.7 ±1 61.2 ± 0.7 ± 1 228 ± 11 ± 24 
46.4 ±1.9 ±2 49.3 ±0.3 ±2 301 ± 3 ± 27 



Table 5 

Summary of the results for the blue data (the first error is the statistical error from 
table 2, and the second the systematic error). 



Epoch 


■^ltt(in m) Aabs (in m) A|^t(in m) 


July 1999 
Sept. 1999 
June 2000 


21.9 ±0.8 ±2 23.5 ±0.1 ±2 119±2±10 
22.8 ±0.3 ±2 25.6 ±0.2 ±2 113±3±10 
26.0 ± 0.5 ±1 28.9 ± 0.1 ±1 133 ± 3 ± 12 



Table 6 

Summary of the results for the UV data (the first error is the statistical error from 
table 3, and the second the systematic error). 

4-6 Stability over the line height 

During the immersion of June 2000, data were recorded with distances of 
100 m and 400 m between the sea bed and the LED source, in both blue and 
UV, to test for a possible variation of the results with depth. As illustrated 
in figure 18, the time distributions are fully compatible with one another, 
whether in blue or in UV, suggesting uniform optical properties along the line 
height. The between the measurements at the two heights is 297 for 330 
degrees of freedom for the blue data, and 252 for 324 degrees of freedom for 
the UV data. 



5 Impact on performance of Antares detector 

The primary goal of the Antares detector is to detect high energy muons 
produced by neutrinos interacting around the detector. In this section, we 
investigate the extent to which the detection of muons in Antares will be 
sensitive to the properties of the water (absorption, scattering and angular 
distribution) and how the uncertainties in the measurement of these properties 
limit the knowledge of the performance of the detector. Since the effect of 
scattering on the time distribution is small, scattering is not expected to have 
a large impact on the reconstruction efficiency, but it might affect the angular 
resolution. 
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Available as a separate figure : comp50_350.jpg 
(a) Time distributions 
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(b) Distributions ratios 



Fig. 18. Comparison of time distributions recorded 100 m above the sea bed (points 
with grey error boxes in (a)) with the distributions recorded 400 m above the sea 
bed (overlaying black line in (a)). Top panels: blue data, bottom panels: UV data. 
The ratio of the distributions, binned into 16 ns bins, are shown in (b). 



5. 1 Event generation and detector simulation 



A sample of 10^^ muon neutrino charged current interactions is generated 
for neutrinos with an E^^''^ spectrum in the energy range 10 GeV < E^ < ?> 
PeV. The angular distribution of the neutrinos is isotropic wdthin the up-going 
hemisphere. The interaction point is uniformly distributed within a cylinder 
of 30 km radius centered on the detector and 30 km height enclosing the 
detector and extending downwards from it. A generator based on LEPTO [23] 
is used for the neutrino interactions. PROPMU [24] is used for those events 
starting outside a 200 m cylinder surrounding the instrumented volume of the 
detector to propagate muons to its surface. Within the 200 m cylinder a full 
detector simulation is then performed including the effect of using different 
scattering models for the photon propagation. The Cherenkov light produced 
by muons and secondary particles is described as a photon field, subsequently 
converted into a photomultiplier hit probability. A final step simulates the 
events in the Antares detector. The detector geometry used is described 
in [25]. The results presented here however are not expected to be strongly 
dependent on the detector configuration and can therefore be applied to the 
present configuration [1]. 

Two parameterisations of the water scattering properties are used. The first 
uses a combination of molecular scattering and tabulated data from "Pet- 



32 



zold" [14] for particle scattering as described in section 3.1. The angular dis- 
tribution is parameterised by a single parameter 77 (cf. equation 5). Two such 
models, PI and P2, are generated with different values of Agct chosen in the 
range of observed values (see table 7), with PI illustrating a conservative 
case. The second parameterisatioiijises a linear combination of two Henyey- 
Greenstein angular distributions (5^^{gi,cos6) to approximate the total scat- 
tering angle distribution: 

;9(cos e)^a cos ^) + (1 - a) /P^(c/2, cos 9) , (18) 

where (3^^ {g, cos 9) is given by equation 16. The parameters describing the 
angular distribution are then a (the relative contribution of the two HG func- 
tions), gl and g2 (the {cos 9) of the two HG functions). Five such models, 
HGl to HG5, are generated with different values of Agct and (cos 6*). 

The simulated models are given in table 7. They have been chosen to include 
both a set of "reasonable" water properties and extreme cases to probe the 
maximum effect of each parameter. The numbers quoted for these models are 
at a wavelength corresponding to the blue LED. The wavelength dependence 
of the scattering lengths is taken according to the Kopelevich model [20]. 



Model 


Asct(m) 


(cos e) 


AL^t(m) 


r] 


a 


gl 


g2 


PI 


40.8 


0.77 


175 


0.17 








P2 


52.0 


0.77 


223 


0.17 








HGl 


52.0 


0.77 


223 




1.000 


0.77 


0.0 


HG2 


22.3 


0.90 


223 




1.000 


0.90 


0.0 


HG3 


4.4 


0.98 


223 




1.000 


0.98 


0.0 


HG4 


40.8 


0.90 


396 




0.985 


0.92 


-0.6 


HG5 


52.0 


0.90 


505 




0.985 


0.92 


-0.6 



Table 7 

Simulated water models and parameters (for A = 466 nm) . 

The scattering angle distributions of all seven models are shown in the upper 
panel of figure 19. The corresponding time distributions 24 m and 44 m away 
from the source, as would be measured with the dedicated ANTARES setup 
(configuration of June 2000 with the blue source), are shown in the lower 
panel. The variety of time distributions is much larger than that actually 
observed in the data (see for comparison figure 14). 

Models PI and P2 include both molecular and particulate contributions to 
the total scattering angle distribution, the molecular part being a major con- 
tributor to the delayed signal, due to its backscattering component which is 
as significant as its forward scattering one. The two models only differ by 
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Available as a separate figure : articleDistrib2.jpg 

Fig. 19. Top panel: angular distributions /?(cos 6) for the models PI, P2, HGl to HG5 
described in the text, as well as for pure molecular and pure particulate scattering 
described in section 3.1 (grey curves). Bottom panel: time distributions for each 
of the models at source-detector distances of 24 m and 44 m (same normalisation 
procedure as in figure 7). From top to bottom (better distinction on the second set 
of curves due to the larger distance of propagation): PI, then P2, then HGl, HG2 
and HG3 almost indistinguishable from one another, then HG4 and finally HG5. 

their scattering lengths. The smaller scattering length of model PI therefore 
generates the larger tail of delayed photons. 

Model HGl reproduces the same Agct and (cos^) as model P2 but with a 
different shape for the scattering angular distribution, as illustrated in the 
upper panel of figure 19. Its lesser backscattering component causes the cor- 
responding time distribution to have less delayed photons. In addition, the 
peak width is increased. Models HG2 and HG3 have the same A^^ as model 
HGl but with different values of (cos^), probing the impact of the angular 
distribution. With HGl, HG2 and HG3 distributions similarly dominated by 
forward scattering, it is the effective scattering length that governs the levels 
of the tail of delayed photons. The corresponding time distributions are thus, 
as expected, almost indistinguishable for the three models. Models HG4 and 
HG5 have the same (cos^) as model HG2 but with different X^^. The angular 
distribution in models HG4 and HG5 is similar to that of model HG2 with 
an enhanced backscattering component. The latter is not sufficient however 
to raise significantly the tail of delayed photons and, as can be expected, the 
increasing effective scattering lengths lowers the levels of the tail of delayed 
photons. 

The absorption profile is the same in all cases and corresponds to that in [3] 
normalised to 62.5 m at 470 nm. 

Model P2 is the one that best reproduces the experimental results described 
in the previous sections. 

5.2 Event reconstruction and analysis 

The three-dimensional reconstruction involves several stages of hit selections 
to remove PMT hits due to "^^K noise and bioluminescence, as well as several 
pre-fits based on plane wave fits through local coincidence and high amplitude 
hits. The final step is based on a maximum likelihood fit to the distribution of 
photon arrival times with respect to the expected arrival time of Cherenkov 
light at a wavelength of 470 nm. The form of the likelihood function is taken 
from an independent Monte Carlo simulation for muons produced by neutri- 
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nos with an E^'^ spectrum above 1 TeV. That Monte Carlo includes a 55 m 
effective attenuation length but neglects scattering so there is no initial as- 
sumption towards a preferred scattering model. Since the likelihood function 
is greatly dominated by the peak where the direct photons were expected, the 
performance of the reconstruction is not strongly affected by the shape of the 
likelihood function. ^° 

The performance of the detector is then defined by two figures of merit of the 
reconstruction: 

Angular resolution Aa defined as the median angle between the Monte 
Carlo neutrino and the reconstructed track. 

Effective volume defined as the fraction of generated events (per bin) which 
remain after reconstruction and selection, multiplied by the generation vol- 
ume. 

Both of these quantities are defined after selection cuts which eliminate mis- 
reconstructed tracks and ensure the purity of the data sample. 

5.3 The effect of different water models 

The angular resolution as defined above is shown for neutrino energies around 
1 TeV (0.3 < <?> TeV) and around 100 TeV (30 < E^ < 300 TeV) for 
each of the simulated water models in figure 20 as a function of the individual 
parameters. Several effects contribute to this resolution. The angle between 
the muon and neutrino at the interaction vertex decreases with increasing 
neutrino energy. At 1 TeV, this angle is 0.7° on average [25] and is the most 
significant contribution to the neutrino angular resolution, whereas at high 
energies the muon and neutrino are essentially collinear so the accuracy of 
reconstructing the muon track dominates the angular resolution. The error 
in the event reconstruction contributes an additional error, bringing the res- 
olution at 1 TeV up to 0.8°. The scattering increases this further to as much 
as 1.2°, depending on the water model. In the high energy regime, relevant 
to neutrino astrophysics, the effect of scattering plays a dominant role. For 
100 TeV, the average muon-neutrino angle is only 0.04°, the reconstruction 
brings it up to 0.20° and the scattering further increases it to as much as 0.53°. 

The remarkable agreement between the angular resolutions obtained for mod- 
els P2 and HGl implies little sensitivity to the precise shape of the scattering 

The event reconstruction and selection described here is not optimal for any 
specific analysis. Rather, it is a general approach with no strong assumptions to 
provide an unbiased assessment of the effect of different scattering models on the 
performance. 
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Angular resolution vs average cos(theta) 



Angular resolution vs scattering length 
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Fig. 20. Angular resolution for each of the water models as a function of 1 — (cos 9), 
Asct and A^^. The upper points (open circles) correspond to the neutrino resolution 
for 0.3 < Ey < ?> TeV with an E~^'^ spectrum, the lower points (filled circles) for 
30 < E,y < 300 TeV. The horizontal lines show the angular resolution obtained in 
the absence of scattering in each case. The effective volumes are shown for the two 
models with the extreme angular resolutions, HG3 (dashed curve) and HG5 (dotted 
curve) and in the absence of scattering (upper solid curve) . 

angular distribution. Same {cos 9) and same Agf^ appear to be sufficient to 
determine the performance of the detector. 

The clearest dependence seen in figure 20 is on the scattering length Asct- A 
shorter scattering length degrades the angular resolution. At 1 TeV, the effect 
of even the most extreme models considered here is at the level of ±12%. 
It is at high energies where the differences between the water models in the 
neutrino angular resolution become most significant, at the level of ±30% 
around the central value. 

The effective scattering length A^f^ alone is not enough to describe the effect 
of scattering on the angular resolution, as virtually the full range of angular 
resolutions obtained are seen for a single effective scattering length of 223 m. 
Similarly, a wide range of values is seen for a single value of {cos 6) = 0.9. 

Given the results of the in-situ measurements (table 2), we can reasonably 
estimate 200 < Af^ < 400 m and a {cos 6) ~ 0.75 for the blue band. With 
this assumption, the variation in angular resolution, even at high energy, is 
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at the level of ±10% around the central value of the angular resolution at 
0.32°. Within this, the effect of assuming different angular distributions (with 
the same {cos 9)) gives an uncertainty at the level of ±6%. Even for the very 
conservative assumption where Agct > 30 m, the uncertainty on the angular 
resolution only increases to 12%. 

The variation in the effective volume of the Antares detector over a wide 
range of neutrino energies between the extreme models is at the level of ±5%, 
with models having the worst angular resolutions also yielding the lowest 
effective volume (see figure 20). 

5.4 Angular resolution of the Antares detector 

The present design of the Antares detector comprises a 12-string network 
that will be immersed over the next few years. The event selection can be 
optimised in terms of angular resolution and effective volume of the detector. 
The track is obtained by a 2-stage fit: the position and orientation of the PMTs 
that have been hit are used to obtain points that the track is likely to have 
crossed. These points are used to obtain an initial track fit. The most probable 
track is then obtained by the minimization of a function involving the residuals 
of the times at which the Cherenkov photons emitted along the track reach the 
PMTs of the detector [26]. In the present stage of the reconstruction software, 
and considering the scattering model that most closely reproduces the data 
presented in this paper, model P2 described above, the angular resolution for 
up-going muon tracks is illustrated in figure 21. For energies > 300 GeV 
the angular resolution for a spectrum is 



The systematics are computed from the study presented in the previous sec- 
tion. 



6 Conclusions 

The fight transmission at the Antares site has been studied intensively with 
dedicated setups designed by the collaboration. Absorption and scattering 

properties of the water for blue light (A = 473 nm) and UV light (A = 375 nm) 
were obtained by measuring the distribution of the arrival times of photons 
emitted by a pulsed LED source and collected several tens of meters away by 
a fast photomultiplier tube. 



Aq;(//) = 0.20° ± 0.01° (stat) ± 0.02° (syst) 
Aa{iy) = 0.32° ± 0.02° (stat) ± 0.04° (syst) . 



(19) 
(20) 
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Fig. 21. Angular resolution (for scattering model P2) as a function of the neutrino 
energy for the reconstruction of the muon track (lower curve) or that of the parent 
neutrino (upper curve) for a 10-string detector. 



The group velocity of light is found in good agreement with predictions from 
analytical models. The absorption length is seen to vary slightly in time, with 
typical values of 60 m in blue and 25 m in UV. These values allow a large 
effective area (~ 0.1 km for neutrino-induced muons with energies in the PeV 
range) for the planned 12-string Antares detector. With the angular distri- 
bution of scattering modelled following the standard approach of oceanogra- 
phers, the scattering length Agct can be extracted with good confidence from 
the data, yielding an effective scattering length A^^ = Asct/(1 — (cos 6*)) of 
~ 260 m in blue and ~ 120 m in UV. The various parameters describing the 
light transmission properties are affected by a 5 to 11% uncertainty, dominated 
by systematics. Given these large scattering lengths, an angular resolution of 
0.3° should be achieved for > 300 GeV, according to the present status 
of the reconstruction software. The uncertainty in the knowledge of the water 
properties (due for instance to the observed variations) affects our knowledge 
of the angular resolution and effective volume of the detector by 10% and 5% 
respectively. 

The light transmission properties were checked to be constant at a given time 
at the two extreme levels (100 m and 400 m above the sea floor) of the active 
part of a detector line. 

The water properties will be monitored with in-situ dedicated instruments 
during the lifetime of the Antares detector so that the instantaneous values 
will be available for use in the muon reconstruction software. 
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